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Abstract 

We include new additional absorption corrections into the Lebiedowicz-Szczurek (non¬ 
resonant) model for pp —>• pprN or pp —)■ pprAn^ processes. They are related to the kN 
nonperturbative interaction in the final state of the reaction. The role of the absorption correc¬ 
tions is quantified for several differential distributions for = 0.2, 1.96, 7, and 8 TeV. The new 
absorption corrections lead to further decrease of the cross section by about a factor of two. They 
change the shape of some distributions {da/dt, da/dpt^p, dcr/dcppp) but leave almost unchanged 
shape of other distributions {da/dMnn, da/dy„, da/dpt^n, da/d(p„j^). The effect may have impor¬ 
tant impact on the interpretation of the recent STAR and CDF data as well as the forthcoming data 
of the ALICE, ATLAS + ALFA and CMS + TOTEM collaborations. 
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I. INTRODUCTION 


There is a growing experimental and theoretical interest in understanding of soft 
hadronic processes at high energy; for reviews see e.g. 10-0] and references therein. 
One of the reaction which can be relatively easy to measure is pp —^ pn^n~p or 
pp —> p7T^TZ~p (four charged particles in the final state). There are recen^ several ex¬ 
perimental projects by the COMPASS 0.01/ STAR 0, CDF ALICE 0[3], ATLAS 
El, and CMS Q collaborations which will measure differential cross sections for the 
reaction(s). Here we wish to compare predictions from the Lebiedowicz-Szczurek model 
with the recent STAR and CDF data. 

The principal reason for studying central exclusive production of mesons is to search 
for glueballs llH . There is some evidence from an analysis of the decay modes of the 
scalar states observed, that the lightest scalar glueball manifests itself through the mixing 
with nearby qq states 11141. IlSll . The exclusive production of lower mass scalar and pseu¬ 
doscalar resonances within a tensor pomeron approach llH was recently examined in 
in . Resonant (p*^ —?► tt+tt”) andnon-resonaritlDrell-Soding) photon-pomeron/reggeon 
production was studied in jl^ . In Refs. 11^0 the continuum background to the pro¬ 
duction of the Tc(O^) state via two-body tt+tt" and K^K~ decays was considered. For 
exclusive production of other mesons see e.g. i2ll-l23l], where mainly the non-central pro¬ 
cesses were discussed. 

Some time ago we proposed a simple phenomenological model for the 7r+7r“- 
continuum mechanism [see diagrams in Fig. [T] (a)] using the tools of Regge theory \2% , 
where perturbative QCD caimot be reliably applied. In this model the parameters of 
pomeron and subleading reggeon exchanges were adjusted to describe total and elastic 
ttN scattering. We expect that our parametrization correctly extrapolate the interaction 
parameters to higher energies where there are no experimental data. The non-resonant 
model can be supplemented to include the pp or pp absorption effects 0 see di¬ 
agrams in Fig. H] (b). The largest uncertainties in the Lebiedowicz-Szczurek model are 
due to the unknown off-shell pion form-factor and the absorption corrections (the soft 
survival factor due to screening corrections). The absorption is done in the eikonal ap¬ 
proximation. The absorption effects lead to substantial damping of the cross section. The 
damping depends on the collision energy and kinematical variables. The discussed here 
model, with reasonable vertex form factors accounting for off-shellness of non-piont-like 
pions in the middle of diagrams in Fig. [T} gives a rough description of the ISR data 


2711 . To get a reasonable description of existing so far experimental data parameter(s) of 
the form factors has (have) to be adjusted 01190. 

It is not clear if the considered so far absorption effects are sufficient to describe the 



FIG. 1: (a) Born amplitudes for the pp —)■ ppn^n process, (b) Absorptive correction amplitudes 
due to the pp interaction. 
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FIG. 2: New absorptive correction amplitudes for the pp —)■ ppn^Ti process due to the zip 
interaction in the final state included in the present analysis. 


data. Any interaction between participating particles potentially leads to absorptive ef¬ 
fects, as it destroys exclusivity of the process. In Fig. [T] (b) we show schematically ab¬ 
sorptive amplitudes due to nN interaction in the final state. Not all combination of in¬ 
teractions are shown in the figure. Some absorptive effects are included inherently in our 
calculation; by using effective interaction fitted to describe nN elastic scattering data. 


Recently, the Lebiedowicz-Szczurek model [24)1 was implemented in GenEx MC |28|| . 
There absorptive corrections are not taken into account explicitly. The authors of [|2^ con¬ 
structed a DIME Monte Carlo code, where almost the same approach was implemented. 
In |l2^ a two-charmel eikonal model was considered. 

In the present paper we will include the additional absorptive effects, not considered 
so far in the literature, and quantify their role for total cross section and tor many differ¬ 
ential distributions for the considered process. Our theoretical results will be compared 
to recent experimental results obtained by the STAR |@] and CDF |@,@] collaborations. We 
will also show some predictions for ALICE, ATLAS and CMS experiments. 


II. BORN AMPLITUDE 

The amplitude squared for the pp —)■ ppn~^n~ process (with four-momenta Pa + Pb ^ 
Pi + P 2 + Ps + P 4 ) considered within the framework of Regge theory with the central 
7r+ n~ system produced by the exchange of two pomeron/reggeon in the t-charmel, as 


shown in Fig. [U, can be written as 

|A4|2 = |A47=oI^ + \Mi=i\^ + \Mi=2\^, (2.1) 

where the isospin amplitudes can be decomposed to the Regge ingredients as 

Mi=o = + {1,0} 1,0 \ 0,0) , (2.2) 

Mi=i = M^P^ + MP^^ + , (2.3) 

M 1=2 = (1,0; 1,0 I 2,0) MP^P^ . (2.4) 

The Clebsch-Gordan coefficients {ji,ini}j 2 ,in 2 \ ],m) are 

(1,0;1,0|0,0) = and (1,0; 1,0 | 2,0) = -VTTs. 

The situation can be summarized as 

|A4j=ol^ ^ |A4/=i|^ 3> |A47=2|^- (2.5) 
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For the dominant pomeron-pomeron contribution we have C-parity C = +1 and isospin 
7 = 0 of the produced 7r~^7r~ system. In general, not only leading double pomeron 
exchanges contribute, but also the subleading f 2 R (C = +1) and pR (C = —1) reggeon 
exchanges. ^ 

The Bom amplitude with the intermediate pion exchange can be written as 

= Mi3(si3, h) M24(s24/ h) + Mu{su, h) M23(S23/ t2)(?..6) 

where the subsystem amplitudes tj) denotes "interaction" between forward pro¬ 

ton (z = 1) or backward proton (z = 2) and one of the two pions {j = 3 for tt"*" or j = 4 
for Tc~). The energy dependence of the np subsystem amplitudes Mij is parametrised in 
terms of the pomeron and the f 2 R reggeon exchange 




rjp Sij 



0Cp{ti)-l 

exp 






exp 



(2.7) 


where z/p = z, = —0.860895 -|- z, Sp is the energy in the (ij) subsystem, and the energy 
scale So is fixed at Sq = 1 GeV^. The pomeron and reggeon trajectories, ap{t) and 
respectively, are assumed to be of standard form, see for instance [1^ , that is, linear in t: 

0 ip{t) = oip{0) +a'pt, ap{0) = 1.0808, a'p = 0.25 GeV“^, (2.8) 

«/„ (0)= 0.5475, = 0.93 GeV-2. (2.9) 

We found the slope parameters fitting the elastic zr+p and n~p cross sections 

= 5.5 GeV“2, = 4 GeV^. (2.10) 

Our model makes various simplifications, but describes the data for elastic nN scattering 
fairly well for energies -s/s^ > 2.5 GeV (see Fig. 2.2 of |0]). 

So far we have assumed a simple exponential dependence of the nN subprocess am¬ 
plitudes (12.7|) which is valid only for small |f| (0.01 < —t < 0.4 GeV^). At larger |t| (ti 
or t 2 in the 2 —> 4 case) the mechanism becomes more complicated. Here a subsequent 
exchange of two pomerons and the exchange of the pomeron together with the reggeon. 


^ The PrPr component is negligible (see the strength parameters in Table 2.1 of ll) and was omitted 
in our analysis. We emphasize, that at lower energies (COMPASS, ISR) the subleading f 2 R exchanges 
constitute a large contribution to the total cross section and must be included in addition to the pomeron 
exchanges; see e.g. section 2.3 of ||3]- Furthermore, there is a large interference effect between the different 
cornponents in the amplitude of about 50% (the total cross section in full phase space), see section 2.6.2 
of As we shall see in the results section imposing limitations on pion rapidity \yTc\ < 1 and going 
to higher energies reduces the role of subleading f 2 R exchanges, however, due to their non-negligible 
interference effects with the leading PP term we keep them explicitly in our calculations. 
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or even pQCD effects (two-gluon exchange) may show up. ^ To give a more realistic t 
dependence we suggest the following replacement (see also 13611 1 


exp 




f{ti,Sij) = exp exp 

B{sij)=Bo + 2a'p\n(^-^^ , 



( 2 . 11 ) 


where the free parameters }i, Bq has been adjusted to the nN elastic scattering data, as 
illustrated in Fig. |3l The so-called 'stretched exponential' parametrization /(t;,s,y) co¬ 
incides at low |f| with the simple exponential form while at larger |f| features a harder 
tail. This function is close to a parametrization ~ exp(—bpj_) suggested by Orear ll^ for 
elastic pp-scattering. The np data show a diffraction dip at — f ~ 4 GeV^. 





FIG. 3: Differential distributions da/dt for zr+p (left panel) and n p (right panel) elastic scattering 
at incident beam momenta Piab = 50 GeV (-^ ~ 9.7 GeV) |38-40] and Pi^b = 200 GeV (-^/s ~ 


19.4 GeV) i38l,l4ll.l42|] . The black dashed lines show results with formula (I2.10D while the red solid 
lines are obtained via the replacement (12.111) . where Bq = 6.5 GeV^^ and = 0.6 GeV^. 


The extra form factors Fr[{t) and Pn{u), in Eq. (12.6)1 . "correct" for the off-shellnes of the 
intermediate pions. The form of the form factor is unknown in particular at higher values 


^ Even description of the elastic pp and pp scattering data is difficult. In Ref. l31ll the authors present a 
model including F + PP + ggg terms and the linear pomeron trajectory. Alternative approaches ll32l. 13311 
combine the soft and hard pomeron exchanges or the odderon exchange in addition 13^ . In the letter 
case the authors consider also a various forms of fhe non-linear pomeron frajecfory. In 13511 fhe role of fhe 
eikonalization of the pp amplitude in both one- and the two-charmel eikonal models were discussed. 


5 































oit or u and they are parametrised in two ways: 


Mt) 

Fnit) 


exp 


Kff,E 

a2 _ ^2 


( 2 . 12 ) 

(2.13) 


and for Fn{u) we have to replace t < — u. These form factors are normalized to unity on 
the pion-mass-shell Fn{in^) = 1. In general, the parameter is not known precisely 
but, in principle, could be fitted to the normalized experimental data. How to extract the 
off-shell parameters will be discussed in the Result section. 


III. ABSORPTION CORRECTIONS 


So far only absorption effects due to pp (pp) interactions were included in the litera¬ 
ture in the calculations of cross sections for the pp —> ppTZ~^7 t~ {pp —)■ ppTz^n~) reactions 
UM- Here we wish to include also the absorption effects due to strong nonperturbative 
interaction of charged pions and (anti)protons in the final state, see corresponding dia¬ 
grams in Fig. |2]). The absorption amplitude including the nN interactions can be written 
in a similar way as that in the case of pp (pp) interaction, i.e., in the eikonal form 


M 


7T p—rescattering 

pp—tpp7T+7T^ 


+ 

+ 

+ 


i 

167r2si4 

i 

167r2si3 

i 

1671^523 

i 

167r2s24 


d% (s, fi, t 2 , ia)M 


P—exchange 

7T~p^7I~p 


(Si4,fc? 


(s, fi, t 2 , ' 


d\t (s, fi, f 2 , h)M 


P—exchange 

7T+p^n+p 


(S23/^? 


j2p Born 

Cl Ki J I pp^ppjj;+ji- 


(S, fi, f2, Wb)A^^-p!J^-"p''(s24, fc?) • (3.1) 


In formula (13.11) we have indicated explicitly only crucial variables, mostly those argu¬ 
ments of Aipp^ppj^+jT:- which get modified in comparison to the Bom amplitude (12.61) . 
For example, the four-momenta squared of the Regge exchange in the first stage of the 
interaction (see Fig. |2]) get modified as 


fl = (Pl - Paf , h = (p2 - Pbf / 


(3.2) 


where the four-momenta of the intermediate nucleons are pi = pi — kt and p 2 = Pi — 
kt- Here, we have introduced auxiliary four-vector kt = (0, kf, 0) to write formulas in a 
compact way. Similarly, the modified four-momenta of pions being propagated in the 
middle of the four-body pp —> ppn^n~ subprocess can be calculated as 

h = {qi - P^f / Ua = {qi - P^f , 

h = {qi - Pif / Wfo = (qi - p^f, (3.3) 


where qi = pa — pi and q2 = Pb — pi are the four-momenta of the (incoming) Regge 
exchanges. We leave all other not explicitly indicated variables which appear in the Born 
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amplitude(s) unchanged. This is an approximation but sufficient for the purpose of the 
present first exploratory analysis. 

The full amplifude includes all rescaffering corrections 


M 


pp—^pp7T~^7T~ 


pp^ppn+n^ 


-U A_^pp-rescattering 
' pp^ppn+ 


-Ur U ^p-rescattering 
' Trp-' pp^ppn^n^ 


(3.4) 


In principle the contributions due to the intermediate proton(s) diffractive excitation(s) 
(p —> N*) could be effectively included by increasing the prefactors. In the present paper 
we shall take, however, Cpp = Cyip = 1. ^ 

In the next section we shall show effect of inclusion of fhe extra absorption terms 
on total cross section as well as on differential distributions. They will lead to further 
decrease of the cross section for the pp —^ ppn~^n~ or pp —> ppn~^TZ~ reactions. We 
expect that the effects may be very important when comparing results of our calculation 
with the recent STAR and CDF experimental data as well as with the forthcoming data of 
the ALICE, CMS, and ATLAS collaborations. 


IV. PREDICTIONS FOR DIFFERENT EXPERIMENTS 

In this section we shall present some selected results for the discussed exclusive pro¬ 
cesses calculated for kinematic domains relevant for the STAR, CDF, ALICE, CMS, and 
ATLAS experiments. In particular, we wish to concentrate on the effect of the new ab¬ 
sorption corrections due to the pion-(anti)proton interaction. We refer also the readers 
to section 2.6.2 of |@] where only the absorptive corrections due to NN interaction were 
discussed. 

Before we go fo fhe higher energies lef us first discuss old ISR data In Fig. H] 

(top panels) we show results for fwo-pion invariant mass distributions. The theoretical 
calculations including absorption corrections have been compared with the ISR data. In 
the calculation the form factor for the off-shell pions was fixed as specified in the fig¬ 
ure captions. The choice of form facfor leads fo differenf behaviour af higher Mnn- We 
also show (bottom panels) the result for fhe exponential and 'stretched exponential' t- 
dependences without and with absorption corrections. The shape of the t distributions 
is strongly modified by the absorption corrections and is similar as obtained in the ISR 
experiment. 


A. STAR experiment 

In Fig. |5] we present the invariant mass distributions of the pion pair produced in 
the pp —> pp7z~^7z~ reaction for fhe STAR kinematics (y^ = 200 GeV, \r]n\ < 1 and 
Pt,n > 0.15 GeV for both pions, the pseudorapidity of the central 7r+7r“ system \r]nn\ < 
2, and in the four-momentum transfers range 0.005 < —fi, —f 2 < 0.03 GeV^). In fhe 
left panel we show result obtained in the Bom approximation (dotted line), the result 
when including proton-proton interactions (dashed line), and when including extra pion- 
nucleon interactions discussed in the present paper. We observe significant damping 


^ How the extra multiplication of the absorption amplitude by a factor Cpp = 1.2 modify 

the features of differential distributions was shown in see, e.g.. Figs. 2.48,2.49,2.50, and Table 2.5. 
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FIG. 4: Two-pion invariant mass distribution at ISR energies with the ISR krnematical cuts indi¬ 
cated in the figure caption. The ISR data I25l.l27|] are shown for comparison. The blue dashed lines 
represent the results obtained for the monopole form factors [ (I2.13L = 1.6 GeV], while the 

black solid lines are for the exponential form r (l2.12L ^ = 1.6 GeV]. The bottom panels repre¬ 

sent the |t| distributions without (dotted lines), with the pp absorption corrections (dashed lines), 
and with all {pp and ttN) absorption corrections included (solid lines). Results for the exponen¬ 
tial (left bottom panel) and for the 'stretched exponential' (right bottom panel) parametrizations 
of the Tip subsystem are shown. 


of the cross section as well as a small shift of the maximum towards smaller invariant 
masses. In the left panel we show result for different parameters of the off-shell form 
factors, together with the STAR experimental data. One can observe that our predictions 
are quite sensitive to the form of the off-shell pion form factor (|2.12|) or (|2.13|) and depend 
on the value of the cut-off parameters If we describe the maximum of the cross 

section around Mnr[ ^ 0.6-0.7 GeV we overestimate the cross section in the interval 1 < 
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Mtttt < 2 GeV essentially for both choices of the form factor form. Part of the effect 
may be related to an enhancement of the cross section due to nn low-energy final state 
interaction. This goes beyond the scope of the present paper which concentrates on the 
new absorption effects. 




M,, (GeV) 


FIG. 5: Two-pion invariant mass distribution at ^/s = 200 GeV with the STAR kinematical cuts 
specified in the figure caption. The dotted line in the left panel corresponds to the Born calcu¬ 
lation, the long-dashed and solid lines to calculations with the absorption effects due to the pp- 
and the 7rp-rescattering in addition, respectively. In the right panel, the blue dashed lines repre¬ 
sent the results with all absorption effects and obtained for the monopole form factors (12.13b for 
different choices of the cut-off parameter Agff^^ = 0.6-1.6 GeV (from bottom to top). The black 
solid lines are for the exponential form (|2.12|l and = 1.0 and 1.6 GeV. The STAR data |0] are 

shown for comparison. 


TT+TT 


In Fig. [6] (left panel) we show differential cross section for the exclusive production of 
system as a function of its pseudorapidity. We conclude that in the range 0.5 < 
Mnn <1-0 GeV both forms of the off-shell pion form factor [ ()2.12|) and (|2.13|) 1 describe 
the data well for = 1.4 — 1.6 GeV. However, the agreement seems a bit misleading 
in the light of disagreement in the invariant mass distribution discussed above. As will 
be discussed in this paper, the absorption effects usually strongly modify the distribution 
in relative azimuthal angle between the outgoing protons (ppp and leave the shape of the 
(pnn distribution almost unchanged. For the STAR (Phase I) visible kinematical range, 
that is at very small four-momentum transfers |f |, one can observe only a damping of the 
cross section (the right bottom panel). The decrease of da/dcppp and da/d(pj 




at (p ^ TT is 


due to the condition \r]j 


< 2 . 


B. CDF experiment 

We wish to emphasize that in this experiment, in contrast to the STAR experiment, 
the final-state nucleons are not detected and only rapidity gap conditions {At] > 4.6 on 
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FIG. 6: The distributions in the pseudorapidity of the produced system and in the 

azimuthal angle between the outgoing pions ((pnn) and between the outgoing protons {(ppp) at 
\/s = 200 GeV in the range of 0.5 < < 1.0 GeV. In the right bottom panel we show the 

results without the absorption effects (the dotted line), with the pp-rescattering (the long-dashed 
line), and with the additional Tip-rescattering (the solid line). The STAR data la] are shown for 
comparison. 


each side of the n~^n~) was imposed experimentally. In Fig. [71 we show the two-pion 
invariant mass distribution at ^/s = 1.96 TeV for the pp —> ppn~^n~ reaction and with 
the following CDF cuts on kinematical variables: pt, 7 T > 0.4 GeV, Ipn] < 1-3 for both 
mesons, and \ymi\ < 1- The rapidity of the central tt+tt" system is expressed by the 
formula 


(P30 + P4o) + (P3z + Plz) \ 

2 V (P30 + P4o) - (P3z + P4z) ) ' 


(4.1) 
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with the four-momenta ps (tt^ meson) and p4 (tt meson). The kinematical cuts pt ,71 ^ 
0.4 GeV on both pions strongly distort the region of low Mnn < 1 GeV. At Mnn — 1 GeV 
the data show a minimum due to interference of the /o(980) resonance contribution with 
the non-resonant background contribution. At higher Mnn, hi the region of 1.2-1.7 GeV, 
some structures could be attributed to /2(1270), /o(1370), /o(1500), and /o(1710) reso¬ 
nant states. The /o(1500) and the /o(1710) mesons are considered to be scalar glueball 
candidates lll3] . but mixing with quarkonium states complicating the issue. We roughly 
describe the differential cross section in the left panel when using the form factors (12.131) 
with Aoff^M — 0-8 GeV. The data at y/s = 0.9 TeV look similar (see Fig. 1 of 1^1]). The data 
at both energies include diffractive dissociation of proton and antiproton (all the pro¬ 
duced unobserved hadrons have | p | > 5.9), so that low diffractive masses of the baryonic 
systems are included, especially at ^/s = 1.96 TeV. In the right panel of Fig. [71 we show 
results with an extra cut on the n~^TZ~ transverse momentum. The results for the form 
factors that give a reasonable description in the left panel badly fail to describe the data 
in the right panel, underestimating the CDF data by a factor of about 5. In this case our 
model results are much below the experimental data which could be due to a contamina¬ 
tion of non-exclusive processes or the perturbative mechanism discussed in |(^ . Both the 
interference of resonant state with the tt"*"/!” -continuum and the diffractive dissociation 
effects require more subtle theoretical approach. This will be addressed elsewhere. Thus, 
the non-resonant Lebiedowicz-Szczurek model should not be expected to fit the data pre¬ 
cisely. We conclude that the CDF data for pt,nn > 1 GeV |@] are sensitive to the details 
of the "large"-f behavior of the nN scattering. While neglecting the extra nN absorptive 
corrections the stretched exponential functional form gives larger cross section than the 
standard exponential form, however, the inclusion of the extra nN absorption cancels the 
"improvement". In this moment we do not know a solution of this puzzle. 

Now let us discussed shortly quantities or observables that are sensitive to the pion 
off-shell form factors. The dependence of < pt,n > and < pt,nn > as a function of two- 
pion invariant mass is presented in Fig. |8l Our calculation shows a rise of the average 
pion transverse momentum with dipion invariant mass. A dependence on the form of 
the form factor is clearly seen. On the contrary, the average transverse momentum of 
the dipion pair is almost independent of the form of the form factor and a parameter of 
the form factor. This can be understood from momentum conservation. The transverse 
momentum of the dipion system must be balanced by the transverse momenta of protons. 
The latter distributions (shapes) are obviously independent of the pion off-shell form 
factors. 

Another observable which can be very sensitive to the choice of off-shell pion form 
factors are the Legendre polynomials < Pieveni^^^^'n+) > (Mnn) distributions, where 

cosO^l' is the angle of the n^ meson with respect to the beam axis, in the 7r+7r“ rest 
frame. In Fig. |9| we present the average Pi calculated as 


< Pl (cos 0^f’) > (Mnn) 


J dVS Pi{cos6^1') da/dVS{Mnn) 
f dVS da/dVS{Mnn) 


(4.2) 


where the integral is done over experimental phase space. We have found that the 
< Pi{cosd^'P) > {Mnn) distributions are almost unaffected by the absorption effects. 
The difference between the results for form factors (|2.12|) and (12.131) is huge at higher 
invariant masses and thus such observables may prove very useful in distinguishing 
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FIG. 7: Two-pion invariant mass distribution at ^/s = 1.96 TeV with the CDF kinematical cuts 
specified in the figure caption. The meaning of the lines in the left panel is the same as in Fig. |5] 
In the right panel we show the results with an additional cut on transverse momentum of the pair 
Pt,nn > 1 GeV without (dotted lines) and with (long-dashed lines) the tzN absorption corrections 
and for two different f-dependences of the Tip-subsystem amplitudes. The black thin lines show 
results with formula (12.10b while the red thick lines represent results for the replacement (12.11)1 
(Bo = 6.5 GeV^^, = 0.6 GeV^). The CDF data |@,S] are shown with only statistical errors; 
systematic uncertainties are approximately 10% at all masses. 



FIG. 8: Mean value of pt,n (left panel) and pt,n 7 z (right panel) as a function of two-pion invariant 
mass at ^/s = 1.96 TeV calculated with the CDF kinematical cuts specified in the figure caption. 

r f 

between these choices. Experimental results of the < Pi{cos6^i') > (Mnn) distri¬ 
butions for ^/s = 1.96 TeV are presented in Fig. 23 of |@] and sfrongly supporf our 
predictions calculafed wifh the monopole form facfors (12.13b and fhe cuf-off paramefer 
Aoff^M = 0-8 GeV, particularly at higher two-pion invariant masses Mnn >1-5 GeV; see 
also discussion in section 2.6.2 (Tevatron) of |S]. One can observe in Fig. 2.42 of 10] thaf 
the contribution of L = 4 is small at low Mnn when the cuts are neglected (left panels) 
and significant already at Mnn ~ 1 GeV when the cuts are applied (right panels). This 
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suggests that the CDF kinematic cuts may distort the partial wave content. This makes 
conclusions more difficult. 



FIG. 9: Mean value of the first even Legendre pol}momials Pi (cos 9 ^(') as a function of two-pion 
invariant mass with the CDF kinematical cuts specified in the figure caption. The results corre¬ 
spond to two types of off-shell pion form factors: the exponential one (12.1211 and the monopole 
one (12.1311 . 


C. ALICE experiment 

Now, we shall present our predictions for experiments at the LFIC. We shall start re¬ 
view of our results for the case of the ALICE experiment at ^/s = 7 TeV. We impose the 
corresponding cuts on both pions transverse momenta pt ,71 ^ 0.1 GeV and pseudora¬ 
pidities \rjTi\ < 0.9. In Fig. [lOlwe show two-pion invariant mass distribution. As for the 
case of the STAR experiment in the left panel we show the Bom result (dotted line), the 
result with pp absorption only (dashed line) as well as the results when including the 
extra zrp absorption (solid line). There is similar tendency as for the STAR case. The ex¬ 
tra absorption lower the cross section without modifying the shape of the invariant mass 
distribution. In the right panel we show our result for two different forms of the off-shell 
form facfor and different values of the cut-off parameters. As for the STAR case the shape 
strongly depends on the form factor form as well as on the cut-off parameters. 

Now we pass to distributions in transverse momenta of single pion and of the pion 
pair, see Fig. [TT] (top panels). The absorption effects due to zrp interaction change the 
shape of the pf^T^Tr-distribution. Such a distribution can be easily measured by the AL¬ 
ICE collaboration. The ALICE experiment caimot register forward/backward protons. 
Therefore only azimuthal correlations between pions can be measured. Our correspond¬ 
ing distribution is shown in FigHT] (bottom panel). The (prcn distribution peaks in the 
back-to-back configuration, i.e. when (prcn = tt. The absorption effects practically do not 
change the shape of the distributions. 

The average values of fransverse momenta of single pion < pt,n > (Mnn) arid of 
the pion pair < pt,nn > (Mnn) are shown in Fig. [121 The results have been obtained 
assuming that pt ,71 ^ 0.1 GeV without the absorption effects (the dotted line), with the 
pp-rescattering (the long-dashed line), and with the additional 7rp-rescattering (the solid 
line). 
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FIG. 10: Two-pion invariant mass distribution at ^/s = 1 TeV with the ALICE kinematical cuts 
specified in the figure caption. In the left panel we show results with a cut-off parameter ^off,E = 
1.6 GeV, see (12.121) . without the absorption effects (the dotted line), with the pp-rescattering (the 
long-dashed line), and with the additional 7rp-rescattering (the solid line). The meaning of the 
lines in the right panel is the same as in Fig. |5l 


D. CMS and ATLAS experiments 

The ATLAS tracking detector provides measurement of charged particle momenta in 
the \r]\ <2.5 region. Since the correlation between the pseudorapidities of both pions is 
very large, the measurement can be performed independently using the tracking detector 
(17/1 < 2.5) and the forward calorimeters (2.5 < \f]\ < 4.9), see Fig. 4 of iHH] . We wish to 
note that the analysis in lUll] was performed for £ = 2 GeV^ neglecting effect of 
the tzN rescattering. Below we shall show results of non-resonant model (including all 
rescattering corrections) for the CMS experiment and the corresponding kinematics cuts 
on both pions: pt^n > 0.1 GeV and \rjr[\ < 2.0. The general features of the differential 
distributions for the ATLAS experiment are, however, similar. 

In Fig. [13] we show two-pion invariant mass distribution. In the left panel we show 
again results for three cases: Bom (dashed line), absorption due to pp interaction (dashed 
line) and for the case with extra np interaction (solid line). In the right panel we show 
the dependence of the cross section on the choice of the pion off-shell form factor. 

Both the CMS (when combined with TOTEM) and ATLAS (when combined with 
ALFA) collaborations can measure outgoing protons. What additional information can 
be provided by measuring the momenta of the outgoing protons? In Figs. [44] and ITSl we 
show the influence of the absorption effects on the p^p, (ppp, and t distributions. The dis¬ 
tribution in proton transverse momenta are particularly interesting. The extra absorption 
effects due to np interactions make the distributions much broader than in the case of 
Born approximation and even broader than in the case when only pp absorption effects 
are included. The effect depends on the value of cut-off parameter Therefore we 

expect that the CMS and ATLAS experimental groups could verify our predictions. The 
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FIG. 11: Differential cross sections dc/d-pt,!!, da/dpt^nn, and dcr/dcp^j,: at ^/s = 1 TeV with the 
ALICE kinematical cuts specified in the figure caption. In the calculation we have used the cut-off 
parameter ^off,E = 1-6 GeV. The meaning of the lines is the same as in Fig. 0 (left panel). 


extra absorption effects lead to significant modification of the shape of proton-proton rel¬ 
ative azimuthal angle distributions which also could be tested by the two experiments. ^ 
The distributions in proton four-momentum transfer t = ti = t 2 are presented in Fig. [14] 
(bottom panels). The extra pion-proton interaction increases the distribution at large |f|. 

The effect of absorption can be even better seen in two-dimensional distributions in 
proton-proton relative azimuthal angle and transverse momentum of one of the protons. 


^ Note that in the case of the ATLAS experiment the requirement of both protons being tagged in the ALFA 
detectors influences the shapes of fhe distributions only very little, but it reduces the cross section by a 
factor close to 3 lU]. 
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FIG. 12: Mean value of the (left panel) and (right panel) as a function of two-pion invari¬ 

ant mass at \/s = 7 TeV with the ALICE kinematical cuts specified in the figure caption. In the 
calculation we have used the cut-off parameter ^off,E = 1-6 GeV. The meaning of the lines is the 
same as in Fig. [TT] 



FIG. 13: Two-pion invariant mass distribution at \/s = 8 TeV with the CMS kinematical cuts 
specified in the figure caption. The meaning of the lines is the same as in Fig. 0 


see [13). Quite different pattern can be seen for the Born case and for the case with full 
absorption. It is not clear to us whether such a two-dimensional distribution can be ob¬ 
tained in practice. 

In Table U we have collected cross sections in fib for the exclusive n^n~ production 
with absorption effects discussed in section |III] and for some kinematical cuts specified 
in section US The Born cross sections for ^/s = 0.2, 1.96, 7, 8 TeV and ^off,E = 1-6 GeV 
are 1.13, 39.60, 54.71, 192.49 fib, respectively. Thus the ratio of full and Bom cross sec¬ 
tions (S^) (the gap survival factor) is approximately 0.20 (STAR), 0.09 (CDF), 0.12 (LHC). 
Results at = 13 TeV were obtained also with the CMS kinematical cuts. 
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TABLE I: The integrated cross sections in jib for the central exclusive 71+production via the 
double-pomeron// 2 R exchange mechanism including the NN and nN absorption effects. The 
results for different experiments with cuts specified in section |lV] and for the different values of 
the off-shell-pion form-factor parameters in Eqs. (12.1211 and (12.1311 are shown. 


Vs (TeV): 

0.2 

1.96 

7 

8 

13 

^off.E = 1-6 GeV 

0.23 

3.69 

6.57 

23.92 

28.64 

^off,E = 1-0 GeV 

0.09 

0.63 

2.16 

7.88 

8.98 

\ff,M = 1-6 GeV 

0.26 

6.45 

9.12 

33.60 

40.92 

^offM = 0.8 GeV 

0.07 

0.58 

1.74 

6.48 

7.45 
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FIG. 14: The distributions in proton transverse momentum (top panels), in azimuthal angle be¬ 
tween the outgoing protons (middle panels), and in proton four-momentum transfer fi (bottom 
panels) at ^/s = 8 TeV with the CMS kinematical cuts specified in the figure caption. In the left 
panels we show the distributions without and with the absorption corrections. In the calculation 
results of which are shown on the right panel we have used two form for the off-shell pion form 
factors and different cut-off parameters Agff. 
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FIG. 15: Two dimensional distributions in pt^p and (ppp at = 8 TeV with the CMS kinematical 
cuts specified in the figure caption. We show the distributions without and with the absorption 
corrections. In the calculation we have used the cut-off parameter ^off,E = 1-6 GeV. 


19 























V. CONCLUSIONS 


In the present paper we have taken into account absorption corrections due to pion- 
nucleon final state interaction in addition to those due to proton-proton interactions. To 
make realistic predictions of the cross sections the parameter responsible for off-shellness 
of intermediate pions in the Lebiedowicz-Szczurek model has been adjusted to experi¬ 
mental data. We have considered here two different scenarios: 

1) The parameters have been adjusted to the STAR data |@], where protons have been 
registered, which guarantees exclusivity of the process, however, the statistics was rather 
low and only low dipion invariant masses (M^zn <1.5 GeV) could be observed. 

2) The parameters have been adjusted to the CDF data |@] (see also |0]), where only some 
rapidity gaps outside of the main detector was imposed in the experiment. 

The cross section for the invariant masses Mjztz < 1 GeV is subjected to low-energy 
pion-pion final state interaction (tttt FSI) effects which are not included in the present 
analysis. Thus, in the first scenario, one finds rather large -^off ~ 1-6 GeV in the region of 
Mjztt < 1 GeV. In the second scenario, when the CDF data are fitted in the broad range of 
Mjztt one obtains ^off ~ 0-8 GeV. Then as a consequence one underestimates the RHIC 
data at M^zn ^ 0.5 — 1.0 GeV. But this missing strerigth at low M^zn is probably due to the 
7T7T FSI enhancement in the cr meson region, see |l43l. 1441] . Therefore, we might expected 
that at higher masses the non-resonant model (with no reggeization of intermediate pion) 
gives realistic predictions with the off-shell pion parameter ~ 1.0 GeV. We have 
proposed to use a 'stretched exponential' parametrization of ttN amplitudes which better 
describes the large-f region and coincides with the exponential parametrization in small- 
t region. Such a parametrization is more adequate when focussing on larger transverse 
momenta. However, we fail to describe the CDF data with pt, 7 Z 7 z > 1 GeV. Clearly final 
tuning of the model requires to take into account both tttt FSI effects as well as explicit 
resonances such as the tensor /2(1270) meson. This goes beyond the scope of the present 
paper, where we have concentrated on the absorption effects. This will be a subject of our 
future studies. 

However, even the present rather simplified treatment of the reaction mechanism al¬ 
lows to draw interesting conclusions as far as the absorption effects are considered. The 
inclusion of the pion-nucleon interactions lead to additional damping of the cross section 
by a factor of about 2, almost independent of center-of-mass energy at least in the range 
considered in the present paper. The additional interaction changes the shape of some 
distributions {da/dti/ 2 , da/dpt,p, da/d(ppp) but leaves almost unchanged shape of other 
distributions {da/dMnn, da/dx/n, da/dpt,n, da/dcpjzn)- Particularly spectacular modifi¬ 
cations are obtained for |f| and pt,p distributions. In particular, a measurement of the 
distribution in the relative azimuthal angle between the pt vectors of the outgoing pro¬ 
tons can provide a fully differential test of the soft survival factors. This could be verified 
in future in experiments when both protons are measured such as ATLAS -t ALFA 111 111 or 
GMS + TOTEM jl^ . In summary, the additional absorption effect discussed here seems 
crucial in detailed understanding of results of ongoing and plaimed experimental inves¬ 
tigations. 
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